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ABSTRACT 

We investigate the rotational equilibrium state of a disk accreting magnetized stars using axisymmetric 
magnetohydrodynamic (MHD) simulations. In this "locked" state, the spin-up torque balances the spin- 
down torque so that the net average torque on the star is zero. We investigated two types of initial 
conditions, one with a relatively weak stellar magnetic field and a high coronal density, and the other 
with a stronger stellar field and a lower coronal density. We observed that for both initial conditions the 
rotation of the star is locked to the rotation of the disk. In the second case, the radial field lines carry 
significant angular momentum out of the star. However, this did not appreciably change the condition 
for locking of the rotation of the star. We find that in the equilibrium state the corotation radius Tco is 
related to the magnetospheric radius ta as Tco/ta ~ 1.2 — 1.3 for case (1) and r coir a « 1.4—1.5 for case 
(2). We estimated periods of rotation in the equilibrium state for classical T Tauri stars, dwarf novae 
and X-ray millisecond pulsars. 

Subject headings: accretion, accretion disks - magnetic fields - plasma - stars: magnetic fields 



1. INTRODUCTION 

Disk accretion to a rotating star with a dipole magnetic 
field is important in a number of astrophysical objects, 
including T Tauri stars (Camenzind 1990; Konigl 1991), 
X-ray pulsars (Nagase 1989; Bildsten et al. 1997), and 
cataclysmic variables (Warner 1995). An important prop- 
erty of this interaction is the disruption of the disk at the 
Alfven radius, ta, and the "locking" of the star's angu- 
lar rotation at an angular velocity, fleq, which is expected 
to be of the order of the disk rotation rate at the Alfven 
radius, fl{rA) — {GM/r^Y^"^. However, the exact con- 
ditions for locking and for the value of the equilibrium 
rotation rate fleq (when the star does not spin up or spin 
down) were not been established. 

One of the complicated aspects of the disk- 
magnetosphere interaction is the process of angular mo- 
mentum transport between the disk and the magnetized 
star. In the first models of the disk-magnetosphere inter- 
action it was proposed that the magnetic field has a dipole 
configuration everywhere and that the net change between 
spin-up torque which arises from the magnetic connection 
of the star to the disk within the corotation radius Vco and 
the spin-down torque which arises from the connection be- 
yond the corotation radius determines the spin evolution 
of the star (Ghosh, Lamb & Pethick 1977; Ghosh & Lamb 
1978, 1979 - hereafter GL79). The corotation radius is the 
radius where the disk rotates at the same speed as a star, 
Tco = {GM/niy/^. It was suggested that for a particular 
value of the star's rotation rate, Qeq, the positive spin-up 
torque balances the negative spin-down torque and the 
star is in the rotational equilibrium state (Ghosh & Lamb 
1978, 1979; Wang 1995). 

Recent studies of the evolution of the magnetic field 
threading the disk and the star led to understanding that 
the field tends to be inflated and possibly opened due to 
the difference in the angular velocities of the foot-points 
(Lovelace, Romanova & Bisnovatyi-Kogan 1995 - hereafter 
LRBK; Shu et al. 1994; Bardou 1999, Uzdensky, Konigl 



& Litwin 2002). In this case a star may lose some angu- 
lar momentum through the open field lines and the equi- 
librium state will be determined by the balance between 
processes of spin- up/spin-down associated with the disk- 
magnetosphere interaction, and spin-down associated with 
the open field lines. 

In some models it was suggested that the angular mo- 
mentum transport between the star and the disk may be 
much less efficient (Agapitou & Papaloizou 2000) if the 
field lines are opened as proposed by LRBK. Under such 
conditions the rotational equilibrium state will be quite 
different (e.g., Matt & Pudritz 2004, 2005). The goal of 
this paper is to derive the conditions for the rotational 
equilibrium state using axisymmetric MHD simulations of 
the disk-magnetosphere interaction. 

The properties of the rotational equilibrium state de- 
pend on the configuration of the magnetic field threading 
the star and the disk. Consequently, analysis of this prob- 
lem requires two or three dimensional simulations. Ax- 
isymmetric simulations have shown that the field lines do 
open (Hayashi et al. 1996; Miller & Stone 1997; Hirose et 
al. 1997; Romanova et al. 1998; Fendt & Elstner 1999). 
However, longer runs have shown that the innermost field 
lines reconnect to form a closed magnetosphere, and some 
of them open and close again in a recurrent manner (Good- 
son k Winglee 1997; Goodson, Winglee & Bohm 1999; 
Romanova et al. 2002, hereafter RUKL; Romanova et 
al. 2004; Kato et al. 2004, Von Rekowski & Branden- 
burg 2004). Detailed simulations of the slow, viscous disk 
accretion to a rotating star with an aligned dipole field 
(RUKL) have shown that on long time-scales, the mag- 
netic field lines in the vicinity of the Alfven radius va are 
closed or only partially open, and these lines are important 
for the angular momentum transport between the star and 
the disk. The balance between the magnetic flux in closed 
and open field lines is clearly important for determining 
the rotational equilibrium state. 

In RUKL, a preliminary search for the conditions of 
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torqueless accretion was performed and the torqueless ac- 
cretion was shown to exist. In this paper, we give a de- 
tailed study of the conditions for torqueless accretion using 
an improved axisymmetric MHD code which makes possi- 
ble longer simulation runs compared to RUKL. The main 
question is: What is the angular rotation rate of a star 
rteq for the torqueless accretion given the star's mass M, 
magnetic moment /i, and accretion rate M . Equivalently, 
if we know the Alfven radius va, then what is the coro- 
tation radius rco in the rotational equilibrium state? In 
earlier theoretical models, it was estimated that the crit- 
ical fastness parameter Ws = r2,/r2ii-(r^) = (rA/fco)^^^ 
of the equilibrium state is in the range of 0.47 — 0.95 (Li 
& Wickramasinghe 1997). This corresponds to the ratio 
rco/rA ~ 1-1 - 1-7. 

In this paper we determine the value of Vco/rA for 
torqueless accretion numerically by performing a large set 
of numerical simulations for different values of 17*, /i, M, 
and disk/corona parameters. We also investigated the 
matter flow and the angular momentum transport in the 
equilibrium state. In §2, the numerical model is described. 
We focus on conditions where the torqueless accretion oc- 
curs and study its dependence on different physical param- 
eters for different cases in §3 and 4. In §5, we apply our 
results to relevant objects, such as CTTSs, cataclysmic 
variables, and X-ray millisecond pulsars. We discuss re- 
sults in Section 6. 



2. THEORETICAL IDEAS AND NUMERICAL SIMULATIONS 



field of the star, ^=B^Rl = Br^ 
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Fig. 1. — Sketch of disk accretion to a star with an ahgned dipole 
magnetic field, where ta is the Alfven radius, rco is corotation ra- 
dius, rip is the deviation radius where a significant deviation from 
Keplerian rotation occurs for r < due to the magnetic force on 
the disk. The dashed line divides the rapidly rotating region and 
the slowly rotating region. FF denotes the funnel flow. 

Theoretical aspects of the disk-magnetoshpere interac- 
tion were developed in the 1970s (e.g. Pringle & Rees, 
1972; Ghosh, Lamb & Pethick 1977; GL79). As sketched 
in Figure 1 , the initial Keplerian accretion disk is threaded 
by the star's magnetic field. The inner radius of the disk, 
referred to as the Alfven radius r^, occurs at the radial 
distance where the energy-density of matter in the disk 
equals to magnetic energy-density, or where the paramter 
/? = pwV(BV87r) = 1. (e.g.. Lamb, Pethick & Pines 1973; 
Davidson & Ostriker 1973). We observed from many simu- 
lations that the inner radius of the disk coincides with this 
radius. The Alfven radius is usually derived analytically 
using a number of approximations. For a dipole magnetic 
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Here, G is the gravitational constant, M is the mass ac- 
cretion rate, M is the mass of the star, r^'' is the Alfven 
radius for the spherical accretion (Eisner & Lamb 1977). 
The coefficient fc^ accounts for the fact that the disk ac- 
cretion is different from spherical accretion. Using values 
of the different parameters obtained from our simulations 
wc obtained r^^^ and then compared it with . From this 
we obtained kA ~ 0.5. 

In a somewhat different approach, Wang (1995) sug- 
gested that the disk stops or disrupts at the radius 
To where the magnetic and matter stresses are equal, 
i-r'^B^Bz)o = M[d{r^Q)/dr]o, where the 0-subscripts 
indicate evaluation at r = tq, B^ and B^ are the compo- 
nents of the magnetic field, fl is the angular velocity of the 
disk. We took values of parameters from our simulations 
and derived the radius the radius tq. We observed that 
this radius approximately coincides with va- 

Whether a star spins up or spins down is due to angu- 
lar momentum flow to or out of the star transported by 
the matter flow and the magnetic field. Some field lines of 
the star thread the disk and transport angular momentum 
between the disk and the star (e.g., GL79). Other field 
lines are infiated and become open, that is, not connected 
to the disk. These open field lines may transport angu- 
lar momentum away from the star to the corona without 
direct interaction with the disk (Agapitou & Papaloizou 
2000; Matt & Pudritz 2004). We observed that both closed 
and open field lines contribute to the angular momentum 
transport to the star. 

For the case where the poloidal field lines connect the 
star and the disk, the field lines passing through the disk 
within the corotation radius rco spin up the star, while 
the field lines passing through the disk at large distances 
spin it down (GL79). The dashed vertical line in Figure 
1 at rco divides the regions which provide positive and 
negative torques. There is also a region r < r$ where the 
magnetic interaction modifies the Keplerian rotation of the 
disk, and where the angular momentum transport between 
the star and the disk is important (GL79, RUKL). If some 
of the magnetic field lines are open, then angular momen- 
tum may be also transported outward from the disk and 
from the star through twisting of open magnetic field lines. 
Thus part of the angular momentum fiux may be trans- 
ported by the matter flow, part by the open field lines, 
and part by closed field lines connecting the star and the 
disk. 

2.1. Main Equations 

The disk-magnetoshpere interaction is considered to be 
described by the MHD equations. 



|+V.(pv) = 0, 



(2) 



P^^-'^P + PS+^i^ xB)xB + F^^\ (3) 
at 47r 



dB 
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V X (v X B) , 
- V • {p^rS) = , 



(4) 
(5) 
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V • B = . (6) 

Here, S = p/ is the entropy per gram, g is the gravita- 
tional acceleration, F"** is the viscous force. The viscosity 
model of Shakura and Sunyaev (1973) is used with viscos- 
ity coefficient v = ac?^/^!^, where c, = {'^p/pY^'^ is the 
sound speed, and parameter a is dimensionless with values 
a~0.01 — 0.03 considered here. 

Equations 2-6 were solved with a Godunov-type numer- 
ical code developed by Koldoba, Kuznetzov & Ustyugova 
(1992) (see also Koldoba & Ustyugova 1994; Ustyugova et 
al. 1995). We used an axisymmetric spherical coordinates 
with a grid x Ng = 131 x 51. Test simulations with 
larger or smaller grids were also performed. 

2.2. Reference Units 

We let Rq denote a reference distance scale, which is 
equal to i?*/0.35 where i?* is the radius of the star. The 
subscript "0" denotes reference values for units. We take 
the reference value for the velocity, vq = {GM/ Rq)^/"^ , the 
angular speed fio = vq/Rq, the timescale Pq = 2-kRo/vq, 
the magnetic field i?o = B^,{R^,/Ro)^ , the magnetic mo- 
ment /xo = BqRq, the density po = Bq/vq, the pressure 
Po = PoVq, the mass accretion rate Mq = PqVoRq, and the 
angular momentum flux Lq = po^o-^o- The dimensionless 
variables are R = R/Rq, v = v/vq, T = T/Pq etc. In the 
following we use these dimensionless units but the tildes 
are implicit. 

2.3. Boundary and Initial Conditions 

Boundary conditions: At the inner boimdary R = R,,, 
"free" boundary conditions are applied for the density, 
pressure, entropy, velocity, and (/)— component of the mag- 
netic field, dp/dR = 0, dp/dR = 0, dS/dR = 0, 
a(v - n X R)i?/9i? = 0, d{RB^)/dR = 0. The poloidal 
components Bj^ and Bg are derived from the magnetic 
flux function ^'(i?, 9), where \1/ at the boundary is derived 
from the frozen-in condition d^/dt + Vp ■ V^* = 0. At the 
outer boundary, free boundary conditions are taken for all 
variables. The outer boundary was placed far away from 
the star in order to diminish the possible influence of this 
boundary (Ustyugova et al. 1999). 

Initial conditions: The star has a fixed aligned dipole 
magnetic field with magnetic moment p. The initial ac- 
cretion disk extends inward to an inner radius r^, and it 
has a temperature Td which is much less than the temper- 
ature of the corona T^. The initial disk and corona are in 
pressure balance. The corona above the disk rotates with 
the angular velocity of the disk in order to avoid a rapid 
initial shearing of the magnetic field. The initial density 
distributions in the disk and corona are constrained by 
the condition that there is a balance of pressure gradient, 
gravitational and centrifuge forces. This balance is neces- 
sary in order to have smooth accretion flow which evolves 
on the viscous time scale of the disk (see details in RUKL) . 
To avoid a very strong initial interaction of the disk with 
the magnetic field, we take the initial inner radius of the 
disk at = 6 or = 5, far away from the star. 

In order to investigate the rotational equilibrium states, 
we performed simulations with two types of initial condi- 
tions, type I and type II. For type I initial conditions we 
have a magnetic field which is relatively weak, jj, = 2, the 



disk is relatively thick with fiducial density pd = 1, and 
relatively high coronal density, poor = 0.005. The initial 
inner disk radius is = 6. The outer radius of simulation 
domain corresponds to 18i?,o, that is, about 54i?*. Results 
for this case are described in §3. For type II initial con- 
ditions the star's magnetic field is much stronger, p — 10, 
the coronal density is much lower pcor = 0.001, the disk 
is thinner, and the initial inner disk radius is = 5. The 
outer radius of simulation domain corresponds to 68i?0i 
that is, about 136i?*. These initial conditions, which are 
more favorable for inflation and opening of the coronal 
fleld lines, are described in §4. 



3. THE EQUILIBRIUM STATE OF DISK ACCRETION FOR 
INITIAL CONDITIONS OF TYPE I 

3.1. Search for Rotational Equilibrium State 

Whether the star spins up or spins down is determined 
by the net flux of angular momentum to the surface of the 
star, L. This flux is composed of two parts, the flux car- 
ried by the matter L^, and that carried by the magnetic 
fleld, Lf. 

L = Lm + Lf, (7) 



Lm = - / dS • prv^Vp , (8) 



Lf = 1- 1 dS- rB^Bp , (9) 

where the p— subscript denotes the poloidal component, 
and dS is the outward pointing surface area element of 
the star. We performed a set of simulations for different 
angular velocities of the star, fi*, to find the critical value 
of n*, corresponding to the rotational equilibrium state, 
that is, the state when L w 0. Other parameters (M, p) 
were fixed. 

We observed that Lf is always dominant over L^ for all 

cases (as in RUKL) and thus we compared Lf for these 
cases. This was predicted in earlier theoretical research 
(e.g., GL79). We narrowed the set of fi* values so that Lb 
was very small on average (see Figure 2a). Among these 
we picked the one for which Lf ^ that corresponds to 
the rotational equilibrium state, fi, — ^eq- For this case, 
the corotation radius {>\.„),;q 1.7. 

In addition, we calculated matter fiux to the star, 

M = - j dS- pvp. (10) 

Figure 2b shows the mentioned fluxes in the rotational 
equilibrium state. Note that this state is typically reached 
at T > 50 because initially the disk is far from the mag- 
netosphere and it takes time for the disk to move inward. 
One can see that the angular momentum flux carried by 
the field Lf fluctuates around zero. 
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Fig. 2. — (a) The evolution of Lf for diflfcrcnt corotation radii 
from Tco = 1-5 (bottom line) to Vco = 2.1 (top line); (b) The evolu- 
tion of M and the angular momentum fluxes carried by matter Lm 
and the magnetic field Lf in rotational equilibrium state, rco ^ 1.7. 
Time scale is Pq (see §2.2). 



3.2. Disk- Magnetosphere Interaction in Equilibrium 

State 

We now discuss in more detail the rotational equilib- 
rium state where rco = 1-7. Figure 3a shows the evolu- 
tion of matter and the magnetic flux with time. One can 
sec that initially the disk matter moves inward, then it 
stops near magnctosphcrc and goes to the star through a 
funnel flow which is driven up by pressure gradient force 
(RUKL;Koldoba et al. 2002). The funnel flow is quasi- 
stationary after a time T ^ 50. The bold red line divides 
the regions where magnetic and matter energy densities 
dominate, that is, the line where (3=1 (see §2). The 
corresponding Alfven radius is ~ 1-2 — 1.3. We also 
estimated the radius ro derived from the equality of mag- 
netic and matter stresses. We observed that this radius is 
close to rA, ro ~ 1-3 — 1.4. 

Figure 3a shows that the magnetic field for r > va 
is strongly non-dipolar, and the structure of the field is 
complicated. The magnetic field lines initially inflate for 
T < 30; however, later some of them reconnect forming 
closed field lines and thus enhancing the dipole component 
in the closed magnetosphere. Other field lines, which are 
near the axis stay open and represent the lines of "mag- 
netic tower" , which is often observed in different simula- 
tions of the disk-magnetosphere interaction (e.g., Kato et 
al. 2004; Romanova et al. 2004). Many field lines are ra- 
dially stretched by the accreting matter. These field lines 
are located in the disk and above the disk. Most of these 



field lines are connected to the star. The field lines above 
the closed magnetosphere continue to open and reconnect 
in quasi-periodic manner. In case of the lower density 
corona (sec §4), inflation is more efficient, and this leads 
to larger quasi-periodic oscillations of the magnetic flux 
and associated fluxes. 

We analyzed the angular momentum transport between 
the star and the disk and corona. Figure 3b shows the 
distribution of the angular momentum fluxes carried by 
the field fs = rB^'Bp/A'K (left panel) and carried by the 
matter = —prv^Vp (right panel). One can see that for 
r > rA, most of the angular momentum flux is carried by 
the matter (see the region with the high positive angu- 
lar momentum at the right panel). However, for r < 
it is mainly transported by the magnetic field (see left 
panel of Figure 3b). The streamlines in the Figure 3b 
show direction of the flow of angular momentum. One can 
see that matter always carries positive angular momentum 
towards the star, which tends to spin up the star. Mag- 
netic fleld lines carry angular momentum out of the star 
through the field lines threading the area of the funnel 
flow and corona. The situation with the angular momen- 
tum transport seems to be more complicated compared to 
one described by GL79. 

The angular momentum fluxes change with distance 
and with angle. We calculated the angular distribu- 
tion of fluxes Fb (r, 6) = 2nr'^ sin 9iB ■ f and Fm (r, 6) =- 
27rr^ sin Q^m'^ through the spheres of different radii r. Fig- 
ure 4 shows that at a large radius, r = 1.8, matter carries 
most of the angular momentum, while the magnetic field 
also contributes but with the opposite sign. At smaller 
radii, the fluxes become smaller and the largest flux is in 
the region of the funnel flow. At the surface of the star, 
the flux associated with matter is very small. The flux as- 
sociated with the field has two components of the opposite 
sign which cancel each other approximately. 

The distribution of angular momentum flux is more 
complicated compared to that of theoretical models. How- 
ever, the equilibrium state does exist and the ratio rco/TA 
is not very much larger than unity. This means that the 
rotation of the star is efficiently locked to the rotation of 
the inner regions of the disk. 

One can see from Figure 3a, that a signiflcant part of 
the disk is disturbed by the disk-magnetosphere interac- 
tion. Figure 5 shows the distribution of the angular veloc- 
ity of the disk in the equatorial plane. We observed that 
the angular velocity varies around the Keplerian value and 
is ustially slightly smaller then Keplerian. The magnetic 
field is strongly wound up in the disk so that the azimuthal 
field dominates. The inner regions of the disk r < ^ 4 
are appreciably influenced by the magnetic fleld. 

3.3. Dependence on fi and a 

Next, wc took a star in the rotational equilibrium state 
(rco = I-''') and varied the magnetic moment of the star 
fi. Figure 6a shows that when ji increases, the rate of 
change of angular momentum of the star becomes more 
negative, that is a star spins-down. This is connected with 
the fact that for larger /j,. the magnetosphere is larger and 
the flux of the positive angular momentum to the star is 
smaller than in the equilibrium state. For example, for 

= 4, the inner disk radius is at « 1.6. At this in- 
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Fig. 3. — (a). Evolution of the density (color background), and the magnetic field (yellow lines) in rotational equilibrium state for type 
I initial conditions for T = — 200 rotations. The density changes from p = 2 in the disk to p = 0.005 in the corona. The bold red line 
corresponds to /3 = 1 (see §2). (b). Fluxes of angular momentum carried by magnetic field (left panel) and by matter fm (right panel) at 
T = 200. Color background shows value of the fiuxcs, the streamlines with arrows - the direction of the angular momentum flow. 



ner radius, the disk rotates much slower than the star and 
a star strongly spins down. Roughly, the dependence is 




Fig. 4. — Angular distribution of fluxes Fm{r,6) and Fg{r,6) for type I initial conditions. The left-hand panel shows radii along which 
the angular momentum fluxes are calculated. The right-hand panel shows the angular momentum fluxes carried by the magnetic fleld (solid 
lines) and the matter (dashed lines) at different radii r. The numbers show the regions where the angular momentum flux carried by the field 
is carried mostly by the closed (I), radial (II), or open (III) field lines. 
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Lf ~ — 0.02/i + 0.04. The balance between spin-up and 
spin-down torques is similar to that suggest by the GL79 
model, where positive and negative torques are associated 
with regions within and beyond the corotation radius. We 
also observed different regions contributing positive and 
negative angular momentum fluxes. However now the dis- 
tribution of the external to the Alfven surface field lines is 
more complicated. 

We calculated the average value of va for each /i and 
obtained the dependence ~ y^*^, with k = 0.36. This 
coefficient k is somewhat different from that of equation 
(1) where k = 4/7 « 0.57. The difference may be con- 
nected with the fact that in the theoretical analysis the 
magnetic field is assumed to be a pure dipole field every- 
where, whereas the simulations show that the actual field 
is different from a dipole for r > r^- In this and all above 
simulations, the viscosity coefficient was fixed at the value 
a = 0.02. 

We also performed simulation runs for different values 
of the disk viscosity, a = 0, 0.01, 0.02, and 0.03, with 
other parameters fixed. We observed that for larger a 
the accretion rate is larger, and the Alfven radius ta is 
smaller. Thus, the region of positive torque becomes larger 
than the region of negative torque, and i/ is larger. In 
other words, at a higher accretion rate M, incoming mat- 
ter brings more of positive angular momentum, which is 
transferred to larger angular momentum flux carried by 
the field Lf. The angular momentum flux increases with 
a (see Figure 6b) as about Lf ~ a^-^ for a > 0.01. 
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Fig. 6. — (a) The evolution of the angular momentum flux to the 
star carried by the magnetic field L f for different magnetic moments 
of the star fi with other parameters fixed, (b) The evolution of the 
angular momentum flux carried by magnetic field Lf for different 
values of of the disk viscosity a with other parameters fixed. 



4. EQUILIBRIUM STATE OF DISK ACCRETION FOR 
INITIAL CONDITIONS OF TYPE II 





Fig. 5. — Radial distribution of angular velocity of the disk along 
in the equatorial plane at T = 180. The dash line shows Keplerian 
angular velocity. 
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T 

Fig. 7. — The evolution of the matter flux M and the angular 
momentum fluxes carried by matter Lm and magnetic fleld i/ in 
the rotational equilibrium state for the type II initial conditions. 

Here we consider another case corresponding to quite 
different initial conditions referred to as type 11. This case 
has a stronger magnetic field, fi — 10, and a much lower 
coronal density, Pcor = 0.001. The density in the corona 
influences the evolution of the magnetic field in the corona, 
namely, the inflation and opening of the magnetic field 
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Fig. 8. — (a). Evolution of the density (color background), and the magnetic field lines (yellow lines) in rotational equilibrium state of type 
II initial conditions for T = — 250 rotations. The density changes from p = 2 in the disk to p = 0.001 in the corona. The bold red line 
corresponds to /? = 1 (see §2). (b). Fluxes of angular momentum carried by magnetic field fg (left panel) and by matter fm (right panel) at 
T = 250. Color background shows value of the fluxes, the streamlines with arrows - the direction of the angular momentum flow. 



lines is favored. As discussed earlier, the angular momen- turn transport between the disk and star occurs in part 




Fig. 9. — Angular distribution of fluxes Fm{r, 9) and Fg{r, 6) for initial conditions of type II. The left-hand panel shows radii along which 
the angular momenta are calculated. The right-hand panel shows angular momentum fluxes carried by magnetic field (solid lines) and matter 
(dashed lines. The numbers show the regions where the angular momentum flux carried by the field is carried mostly by the closed (I), radial 
(II), or open (III) field lines. 
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through the magnetic field. In addition, we took a dif- 
ferent structure of the disk, which is thinner and located 
initially at = 5, which is closer to the star for such a 
strong field. Wc were able to model the stronger magnetic 
field because wc increased the radius of the star (the in- 
ner boundary) to i?* = 0.5 (versus i?, = 0.35 in other 
simulation runs). In addition, we were able to use a grid 
Nr X Ng = 91x31 which speeded up the simulations and 
allowed longer nms, up to 1000 Pq. 

For this case, we also performed a set of simulations for 
different corotation radii and found the state which gives 
torqueless accretion, L « 0. We found that this state cor- 
responds to a corotation radius Tco ~ 4.6. Figure 7 shows 
the evolution of the fluxes in this case. The angular mo- 
mentum flux associated with the matter is very small as 
in our main case. The flux due to the magnetic field, Lf, 
varies more strongly than in the main case. This is be- 
cause the coronal density is lower so that variations of the 
magnetic field in the corona are larger. In addition, the 
accretion rate M is several times larger than that for our 
main case, because the disk is different compared to the 
main case. 

We analyzed the rotational equilibrium state in detail. 
Figure 8a shows the evolution of the disk and magnetic 
field in the inner part of the simulation region. We ob- 
served that the accretion disk stopped at much larger 
distances, ta ~ 3.1, compared to our main case (r^ ~ 
1.2 — 1.3) because of much stronger magnetic field. This 
radius again coincides with the (3 = 1 line (see the bold 
red line on figure 8a). Wc also estimated the radius ro 
derived from comparison of magnetic and matter torques. 
Again, we obtained similar but a somewhat larger radius, 
ro « 3.2. We observed that in type II case, the ratio 
fco/TA ^ 1-4 — 1.5. This number is slightly larger than 
that in the main case {tco/ta ^ 1-3 — 1.4) which may be 
connected with contribution of the inflated field lines to 
the negative torque. 

Figure 8b shows the angular momentum fluxes carried 
by the field -Fs(r, 9) and by the matter Fm(r, 9) . One can 
see that the disk matter brings positive spin-up torque 
(right panel), while magnetic field carries both negative 
and positive torques (left panel). In this case a significant 
part of the spin-down fiux is carried by the inflated and ra- 
dially stretched field lines. This is different from the type 
I initial conditions case. 

Figure 9 shows angular distribution of ffuxes Fb {r, 6) 
and Fm{r,9) along the spheres of different radii r from 
r = 0.5 to r = 4.5. One can see, that at large distances, 
r = 4.5, there is a large positive angular momentum flux 
carried by the matter of the disk, and there is a large 
negative angular momentum flux carried by the radially 
stretched field lines located above the disk and in the disk. 
Field lines above the disk are inflated field lines which con- 
nect to the disk at large radii r > 20. They carry both, 
positive matter flux, and negative spin-down flux. In this 
region the density is smaller than in the disk, however, 
the field lines are strongly wound and azimuthal velocity 
is large, this is why both fluxes associated with matter 
and the field are large in this region. The plot at r = 3.75 
shows relative input from different sets of field lines. It 
shows that in the region of the closed field lines (region I) 
matter carries positive torque while field lines carry only 



small negative torque. Most of both torque is carried by 
the radially stretched field lines (region II), while the role 
of vertically inflated field lines of "magnetic tower" (region 
III) is very small. At even smaller radii, r = 3,2,1 and 
0.5, both torques become smaller, like in type I case. 

5. EQUILIBRIUM STATE IN APPLICATIONS TO DIFFERENT 

STARS 

The rotational equilibrium state is expected to be the 

most probable state for different disk accreting magnetized 
stars. These include Classical T Tauri Stars, cataclysimic 
variables, and X-ray pulsars. For example, the slowly ro- 
tating CTTSs have spin rates as low as ~ 10% of breakup 
speed (Bouvier et al. 1993). The fact that CTTSs ro- 
tate slowly strongly suggests that they are in rotational 
equilibrium state. 

Below we estimate periods of rotation for different 
weakly magnetized stars. Our simulations were performed 
for cases ta/R* =4 — 6 and can be applied to systems 
with different scales in which this ratio is satisfied. Ob- 
servations show that in CTTS this ratio is in the range 
r^/i?^ = 4-8 (Bouvier et. al., 1993, Edwards et al., 1993), 
or, rA/R* = 3— 10 (Kenyon, Yi & Hartmann, 1996). Sim- 
ilar or smaller magnetospheres are expected in accreting 
millisecond pulsars (van der Klis 2000) and dwarf novae 
cataclysmic variables (Warner 2004). 

For our example we take va/R* = 4, take obtained from 
simulation ratio rco = I-'Ita and derive the angular veloc- 
ity for rotational equilibrium, Oe^: 

3/2 

w 0.09 , (11) 
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which is about ~ 10% of breakup speed of the star. Here, 
il^K is the Keplerian angular speed at i?*. The value of 
*if we find is close to the one observed in CTTSs. 

5.1. Classical T Tauri Stars (CTTSs) 

The angular velocity of the star can be written as 
ri* = {GM /r^gY^"^ . Combining this with the rotational 
equilibrium condition rco/rA ^ 1-3 — 1.5 (an average 1.4 
was taken here) and using Eqn. 1 for the Alfven radius, 
we can get the critical value of angular velocity Ogg or ro- 
tation period Peq = 2n/^eq- For the CTTS we obtain the 
equilibrium rotation period, 

P.,.4.6d.ys(!^)"'(i2::|^)="x 

where B is the surface magnetic field of the star. Our re- 
sult is in the range of an observational bimodal distribution 
of rotation at period with peaks near 3 and 8 days (At- 
tridge & Herbst 1992). Edwards et al. (1993) proposed 
that the distribution near 8 days may be caused by the 
disk locking. We should note that both peaks may be con- 
nected with the disk locking, but at different parameters, 
say, M and B. 

5.2. Cataclysmic Variables 

We next consider dwarf novae which belong to a sub- 
type of the cataclysmic variables where the magnetic field 
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is expected to be small but still possibly dynamically im- 
portant for the accretion. The accretion disks are stopped 
at small radii by the white dwarf magnetosphere. The ac- 
creting material then leaves the disk and follows the mag- 
netic field lines down to the star's surface in the vicinity 
of the magnetic poles. Taking typical values for the white 
dwarf and accretion disk, we obtain the period of these 
stars in rotational equilibrium state, 



eg 



(■ 



M 
R 



X 18/7 

VlQSG/ V7.0xl08cm7 ' ^^^^ 

The observed periods of Dwarf Novae Oscillations (DNOs) 
are in the range of 7 — 70 s (Warner, 2004). 

5.3. Millisecond Pulsars 

Our simulations are also applicable to millisecond X-ray 
pulsars, which are the accreting low-magnetic-field neu- 
tron stars (van dcr Klis, 2000). The period of such a star 
in the rotational equilibrium is 

-(^)"'(i^)"'^ 

/ B \6/7/ R y8/7 

.lO^cm/ 



IiOSqJ Uo6cmJ ■ ^^^^ 
For all accretion-powered millisecond X-ray pulsars we 

have known, the spin frequencies range from 185Hz to 
435Hz (Wijnands, R. 2004), or, the spin rate is in range 
of 2.3-5.2 ms. These values may be obtained from Eqn.l4 
at smallere M or larger B. 

6. DISCUSSION 

We investigated the conditions for the rotational equi- 
librium or"torqueless" accretion state using axisymmetric 
MHD simulations. In such a state the total angular mo- 
mentum flux to the star is approximately zero. We found 
the equilibrium states by gradually changing the angu- 
lar velocity of the star with the other parameters fixed. 
We considered two main cases: One with relatively low 
magnetic field and dense corona and second with much 
stronger field and lower density corona. We observed that 
in both cases the rotation of the star is approximately 
locked to the rotation of the inner radius of the disk such 



that a star rotates somewhat slower than the inner radius 
of the disk. In the first case the ratio between corota- 
tion radius of the star and the Alfven radius (where 
the disk is disrupted) is rco/rA ^ 1-3 — 1.4. In the second 
case, where stronger inflation of the magnetic field was ob- 
served, tliis ratio is only slightly larger, rco/rA ^ 1.4— 1.5. 
We observed that in the first case, the angular momen- 
tum transport is associated with the closed field lines at 
the inner radii of the disk. Open field lines spin down the 
star, but the role of this spin-down is small. In the sec- 
ond case outflow of angular momentum along the inflated 
field lines is more significant than in the first case, how- 
ever this did not change the result, which probably means 
that the angular momentum transport associated with the 
inner regions of the disk and the region of the funnel flow 
is more significant. Thus, in both cases the magnetic in- 
teraction effectively locks the rotation rate of the star to 
a value which depends mainly on the mass accretion rate 
and the star's magnetic moment. We should note that 
most of coronal region is still matter dominated. This is 
connected with the fact that the magnetic energy-density 
of the dipole decreases with distance as ^ R^^ so that 
it is difficult to set up a magnetically dominated corona. 
Goodson et al. (1997, 1999) were able to model such a 
corona by arranging fast fall-off of coronal density with 
the distance. However, their initial conditions are sufS- 
cicntly far from equilibrium that the torqueless accretion 
was not established. Future simulations with even lower 
coronal densities will help to understand whether a star 
is always locked to the rotation rate of the inner radius of 
the disk. Our simulations for two very different initial con- 
ditions have shown very similar results for the rotational 
equilibrium which may be a sign that the disk-locking may 
be a similar for all slowly rotating stars. 

We applied our simulation results to Classical T Tauri 
stars, where disk locking may explain their slow rotation. 
Also, we estimated the probable periods of rotation of 
other accretion powered systems, such as dwarf novae and 
X-ray millisecond pulsars. 
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ABSTRACT 

We investigate the rotational equilibrium state of a disk accreting magnetized stars using axisymmetric 
magnetohydrodynamic (MHD) simulations. In this "locked" state, the spin-up torque balances the spin- 
down torque so that the net average torque on the star is zero. We investigated two types of initial 
conditions, one with a relatively weak stellar magnetic field and a high coronal density, and the other 
with a stronger stellar field and a lower coronal density. We observed that for both initial conditions the 
rotation of the star is locked to the rotation of the disk. In the second case, the radial field lines carry 
significant angular momentum out of the star. However, this did not appreciably change the condition 
for locking of the rotation of the star. We find that in the equilibrium state the corotation radius Tco is 
related to the magnetospheric radius ta as Tco/ta ~ 1.2 — 1.3 for case (1) and r coir a « 1.4—1.5 for case 
(2). We estimated periods of rotation in the equilibrium state for classical T Tauri stars, dwarf novae 
and X-ray millisecond pulsars. 

Subject headings: accretion, accretion disks - magnetic fields - plasma - stars: magnetic fields 



1. INTRODUCTION 

Disk accretion to a rotating star with a dipole magnetic 
field is important in a number of astrophysical objects, 
including T Tauri stars (Camenzind 1990; Konigl 1991), 
X-ray pulsars (Nagase 1989; Bildsten et al. 1997), and 
cataclysmic variables (Warner 1995). An important prop- 
erty of this interaction is the disruption of the disk at the 
Alfven radius, ta, and the "locking" of the star's angu- 
lar rotation at an angular velocity, fleq, which is expected 
to be of the order of the disk rotation rate at the Alfven 
radius, fl{rA) — {GM/r^Y^"^. However, the exact con- 
ditions for locking and for the value of the equilibrium 
rotation rate fleq (when the star does not spin up or spin 
down) were not been established. 

One of the complicated aspects of the disk- 
magnetosphere interaction is the process of angular mo- 
mentum transport between the disk and the magnetized 
star. In the first models of the disk-magnetosphere inter- 
action it was proposed that the magnetic field has a dipole 
configuration everywhere and that the net change between 
spin-up torque which arises from the magnetic connection 
of the star to the disk within the corotation radius Vco and 
the spin-down torque which arises from the connection be- 
yond the corotation radius determines the spin evolution 
of the star (Ghosh, Lamb & Pethick 1977; Ghosh & Lamb 
1978, 1979 - hereafter GL79). The corotation radius is the 
radius where the disk rotates at the same speed as a star, 
Tco = {GM/niy/^. It was suggested that for a particular 
value of the star's rotation rate, Qeq, the positive spin-up 
torque balances the negative spin-down torque and the 
star is in the rotational equilibrium state (Ghosh & Lamb 
1978, 1979; Wang 1995). 

Recent studies of the evolution of the magnetic field 
threading the disk and the star led to understanding that 
the field tends to be inflated and possibly opened due to 
the difference in the angular velocities of the foot-points 
(Lovelace, Romanova & Bisnovatyi-Kogan 1995 - hereafter 
LRBK; Shu et al. 1994; Bardou 1999, Uzdensky, Konigl 



& Litwin 2002). In this case a star may lose some angu- 
lar momentum through the open field lines and the equi- 
librium state will be determined by the balance between 
processes of spin- up/spin-down associated with the disk- 
magnetosphere interaction, and spin-down associated with 
the open field lines. 

In some models it was suggested that the angular mo- 
mentum transport between the star and the disk may be 
much less efficient (Agapitou & Papaloizou 2000) if the 
field lines are opened as proposed by LRBK. Under such 
conditions the rotational equilibrium state will be quite 
different (e.g., Matt & Pudritz 2004, 2005). The goal of 
this paper is to derive the conditions for the rotational 
equilibrium state using axisymmetric MHD simulations of 
the disk-magnetosphere interaction. 

The properties of the rotational equilibrium state de- 
pend on the configuration of the magnetic field threading 
the star and the disk. Consequently, analysis of this prob- 
lem requires two or three dimensional simulations. Ax- 
isymmetric simulations have shown that the field lines do 
open (Hayashi et al. 1996; Miller & Stone 1997; Hirose et 
al. 1997; Romanova et al. 1998; Fendt & Elstner 1999). 
However, longer runs have shown that the innermost field 
lines reconnect to form a closed magnetosphere, and some 
of them open and close again in a recurrent manner (Good- 
son k Winglee 1997; Goodson, Winglee & Bohm 1999; 
Romanova et al. 2002, hereafter RUKL; Romanova et 
al. 2004; Kato et al. 2004, Von Rekowski & Branden- 
burg 2004). Detailed simulations of the slow, viscous disk 
accretion to a rotating star with an aligned dipole field 
(RUKL) have shown that on long time-scales, the mag- 
netic field lines in the vicinity of the Alfven radius va are 
closed or only partially open, and these lines are important 
for the angular momentum transport between the star and 
the disk. The balance between the magnetic flux in closed 
and open field lines is clearly important for determining 
the rotational equilibrium state. 

In RUKL, a preliminary search for the conditions of 
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torqueless accretion was performed and the torqueless ac- 
cretion was shown to exist. In this paper, we give a de- 
tailed study of the conditions for torqueless accretion using 
an improved axisymmetric MHD code which makes possi- 
ble longer simulation runs compared to RUKL. The main 
question is: What is the angular rotation rate of a star 
rteq for the torqueless accretion given the star's mass M, 
magnetic moment /i, and accretion rate M . Equivalently, 
if we know the Alfven radius va, then what is the coro- 
tation radius rco in the rotational equilibrium state? In 
earlier theoretical models, it was estimated that the crit- 
ical fastness parameter Ws = r2,/r2ii-(r^) = (rA/fco)^^^ 
of the equilibrium state is in the range of 0.47 — 0.95 (Li 
& Wickramasinghe 1997). This corresponds to the ratio 
rco/rA ~ 1-1 - 1-7. 

In this paper we determine the value of Vco/rA for 
torqueless accretion numerically by performing a large set 
of numerical simulations for different values of 17*, /i, M, 
and disk/corona parameters. We also investigated the 
matter flow and the angular momentum transport in the 
equilibrium state. In §2, the numerical model is described. 
We focus on conditions where the torqueless accretion oc- 
curs and study its dependence on different physical param- 
eters for different cases in §3 and 4. In §5, we apply our 
results to relevant objects, such as CTTSs, cataclysmic 
variables, and X-ray millisecond pulsars. We discuss re- 
sults in Section 6. 



2. THEORETICAL IDEAS AND NUMERICAL SIMULATIONS 



field of the star, ^=B^Rl = Br^ 
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Fig. 1. — Sketch of disk accretion to a star with an ahgned dipole 
magnetic field, where ta is the Alfven radius, rco is corotation ra- 
dius, rip is the deviation radius where a significant deviation from 
Keplerian rotation occurs for r < due to the magnetic force on 
the disk. The dashed line divides the rapidly rotating region and 
the slowly rotating region. FF denotes the funnel flow. 

Theoretical aspects of the disk-magnetoshpere interac- 
tion were developed in the 1970s (e.g. Pringle & Rees, 
1972; Ghosh, Lamb & Pethick 1977; GL79). As sketched 
in Figure 1 , the initial Keplerian accretion disk is threaded 
by the star's magnetic field. The inner radius of the disk, 
referred to as the Alfven radius r^, occurs at the radial 
distance where the energy-density of matter in the disk 
equals to magnetic energy-density, or where the paramter 
/? = pwV(BV87r) = 1. (e.g.. Lamb, Pethick & Pines 1973; 
Davidson & Ostriker 1973). We observed from many simu- 
lations that the inner radius of the disk coincides with this 
radius. The Alfven radius is usually derived analytically 
using a number of approximations. For a dipole magnetic 
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Here, G is the gravitational constant, M is the mass ac- 
cretion rate, M is the mass of the star, r^'' is the Alfven 
radius for the spherical accretion (Eisner & Lamb 1977). 
The coefficient fc^ accounts for the fact that the disk ac- 
cretion is different from spherical accretion. Using values 
of the different parameters obtained from our simulations 
wc obtained r^^^ and then compared it with . From this 
we obtained kA ~ 0.5. 

In a somewhat different approach, Wang (1995) sug- 
gested that the disk stops or disrupts at the radius 
To where the magnetic and matter stresses are equal, 
i-r'^B^Bz)o = M[d{r^Q)/dr]o, where the 0-subscripts 
indicate evaluation at r = tq, B^ and B^ are the compo- 
nents of the magnetic field, fl is the angular velocity of the 
disk. We took values of parameters from our simulations 
and derived the radius the radius tq. We observed that 
this radius approximately coincides with va- 

Whether a star spins up or spins down is due to angu- 
lar momentum flow to or out of the star transported by 
the matter flow and the magnetic field. Some field lines of 
the star thread the disk and transport angular momentum 
between the disk and the star (e.g., GL79). Other field 
lines are infiated and become open, that is, not connected 
to the disk. These open field lines may transport angu- 
lar momentum away from the star to the corona without 
direct interaction with the disk (Agapitou & Papaloizou 
2000; Matt & Pudritz 2004). We observed that both closed 
and open field lines contribute to the angular momentum 
transport to the star. 

For the case where the poloidal field lines connect the 
star and the disk, the field lines passing through the disk 
within the corotation radius rco spin up the star, while 
the field lines passing through the disk at large distances 
spin it down (GL79). The dashed vertical line in Figure 
1 at rco divides the regions which provide positive and 
negative torques. There is also a region r < r$ where the 
magnetic interaction modifies the Keplerian rotation of the 
disk, and where the angular momentum transport between 
the star and the disk is important (GL79, RUKL). If some 
of the magnetic field lines are open, then angular momen- 
tum may be also transported outward from the disk and 
from the star through twisting of open magnetic field lines. 
Thus part of the angular momentum fiux may be trans- 
ported by the matter flow, part by the open field lines, 
and part by closed field lines connecting the star and the 
disk. 

2.1. Main Equations 

The disk-magnetoshpere interaction is considered to be 
described by the MHD equations. 



|+V.(pv) = 0, 



(2) 



P^^-'^P + PS+^i^ xB)xB + F^^\ (3) 
at 47r 



dB 

'dt ' 
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V • B = . (6) 

Here, S = p/ is the entropy per gram, g is the gravita- 
tional acceleration, F"** is the viscous force. The viscosity 
model of Shakura and Sunyaev (1973) is used with viscos- 
ity coefficient v = ac?^/^!^, where c, = {'^p/pY^'^ is the 
sound speed, and parameter a is dimensionless with values 
a~0.01 — 0.03 considered here. 

Equations 2-6 were solved with a Godunov-type numer- 
ical code developed by Koldoba, Kuznetzov & Ustyugova 
(1992) (see also Koldoba & Ustyugova 1994; Ustyugova et 
al. 1995). We used an axisymmetric spherical coordinates 
with a grid x Ng = 131 x 51. Test simulations with 
larger or smaller grids were also performed. 

2.2. Reference Units 

We let Rq denote a reference distance scale, which is 
equal to i?*/0.35 where i?* is the radius of the star. The 
subscript "0" denotes reference values for units. We take 
the reference value for the velocity, vq = {GM/ Rq)^/"^ , the 
angular speed fio = vq/Rq, the timescale Pq = 2-kRo/vq, 
the magnetic field i?o = B^,{R^,/Ro)^ , the magnetic mo- 
ment /xo = BqRq, the density po = Bq/vq, the pressure 
Po = PoVq, the mass accretion rate Mq = PqVoRq, and the 
angular momentum flux Lq = po^o-^o- The dimensionless 
variables are R = R/Rq, v = v/vq, T = T/Pq etc. In the 
following we use these dimensionless units but the tildes 
are implicit. 

2.3. Boundary and Initial Conditions 

Boundary conditions: At the inner boimdary R = R,,, 
"free" boundary conditions are applied for the density, 
pressure, entropy, velocity, and (/)— component of the mag- 
netic field, dp/dR = 0, dp/dR = 0, dS/dR = 0, 
a(v - n X R)i?/9i? = 0, d{RB^)/dR = 0. The poloidal 
components Bj^ and Bg are derived from the magnetic 
flux function ^'(i?, 9), where \1/ at the boundary is derived 
from the frozen-in condition d^/dt + Vp ■ V^* = 0. At the 
outer boundary, free boundary conditions are taken for all 
variables. The outer boundary was placed far away from 
the star in order to diminish the possible influence of this 
boundary (Ustyugova et al. 1999). 

Initial conditions: The star has a fixed aligned dipole 
magnetic field with magnetic moment p. The initial ac- 
cretion disk extends inward to an inner radius r^, and it 
has a temperature Td which is much less than the temper- 
ature of the corona T^. The initial disk and corona are in 
pressure balance. The corona above the disk rotates with 
the angular velocity of the disk in order to avoid a rapid 
initial shearing of the magnetic field. The initial density 
distributions in the disk and corona are constrained by 
the condition that there is a balance of pressure gradient, 
gravitational and centrifuge forces. This balance is neces- 
sary in order to have smooth accretion flow which evolves 
on the viscous time scale of the disk (see details in RUKL) . 
To avoid a very strong initial interaction of the disk with 
the magnetic field, we take the initial inner radius of the 
disk at = 6 or = 5, far away from the star. 

In order to investigate the rotational equilibrium states, 
we performed simulations with two types of initial condi- 
tions, type I and type II. For type I initial conditions we 
have a magnetic field which is relatively weak, jj, = 2, the 



disk is relatively thick with fiducial density pd = 1, and 
relatively high coronal density, poor = 0.005. The initial 
inner disk radius is = 6. The outer radius of simulation 
domain corresponds to 18i?,o, that is, about 54i?*. Results 
for this case are described in §3. For type II initial con- 
ditions the star's magnetic field is much stronger, p — 10, 
the coronal density is much lower pcor = 0.001, the disk 
is thinner, and the initial inner disk radius is = 5. The 
outer radius of simulation domain corresponds to 68i?0i 
that is, about 136i?*. These initial conditions, which are 
more favorable for inflation and opening of the coronal 
fleld lines, are described in §4. 



3. THE EQUILIBRIUM STATE OF DISK ACCRETION FOR 
INITIAL CONDITIONS OF TYPE I 

3.1. Search for Rotational Equilibrium State 

Whether the star spins up or spins down is determined 
by the net flux of angular momentum to the surface of the 
star, L. This flux is composed of two parts, the flux car- 
ried by the matter L^, and that carried by the magnetic 
fleld, Lf. 

L = Lm + Lf, (7) 



Lm = - / dS • prv^Vp , (8) 



Lf = 1- 1 dS- rB^Bp , (9) 

where the p— subscript denotes the poloidal component, 
and dS is the outward pointing surface area element of 
the star. We performed a set of simulations for different 
angular velocities of the star, fi*, to find the critical value 
of n*, corresponding to the rotational equilibrium state, 
that is, the state when L w 0. Other parameters (M, p) 
were fixed. 

We observed that Lf is always dominant over L^ for all 

cases (as in RUKL) and thus we compared Lf for these 
cases. This was predicted in earlier theoretical research 
(e.g., GL79). We narrowed the set of fi* values so that Lb 
was very small on average (see Figure 2a). Among these 
we picked the one for which Lf ^ that corresponds to 
the rotational equilibrium state, fi, — ^eq- For this case, 
the corotation radius {>\.„),;q 1.7. 

In addition, we calculated matter fiux to the star, 

M = - j dS- pvp. (10) 

Figure 2b shows the mentioned fluxes in the rotational 
equilibrium state. Note that this state is typically reached 
at T > 50 because initially the disk is far from the mag- 
netosphere and it takes time for the disk to move inward. 
One can see that the angular momentum flux carried by 
the field Lf fluctuates around zero. 
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Fig. 2. — (a) The evolution of Lf for diflfcrcnt corotation radii 
from Tco = 1-5 (bottom line) to Vco = 2.1 (top line); (b) The evolu- 
tion of M and the angular momentum fluxes carried by matter Lm 
and the magnetic field Lf in rotational equilibrium state, rco ^ 1.7. 
Time scale is Pq (see §2.2). 



3.2. Disk- Magnetosphere Interaction in Equilibrium 

State 

We now discuss in more detail the rotational equilib- 
rium state where rco = 1-7. Figure 3a shows the evolu- 
tion of matter and the magnetic flux with time. One can 
sec that initially the disk matter moves inward, then it 
stops near magnctosphcrc and goes to the star through a 
funnel flow which is driven up by pressure gradient force 
(RUKL;Koldoba et al. 2002). The funnel flow is quasi- 
stationary after a time T ^ 50. The bold red line divides 
the regions where magnetic and matter energy densities 
dominate, that is, the line where (3=1 (see §2). The 
corresponding Alfven radius is ~ 1-2 — 1.3. We also 
estimated the radius ro derived from the equality of mag- 
netic and matter stresses. We observed that this radius is 
close to rA, ro ~ 1-3 — 1.4. 

Figure 3a shows that the magnetic field for r > va 
is strongly non-dipolar, and the structure of the field is 
complicated. The magnetic field lines initially inflate for 
T < 30; however, later some of them reconnect forming 
closed field lines and thus enhancing the dipole component 
in the closed magnetosphere. Other field lines, which are 
near the axis stay open and represent the lines of "mag- 
netic tower" , which is often observed in different simula- 
tions of the disk-magnetosphere interaction (e.g., Kato et 
al. 2004; Romanova et al. 2004). Many field lines are ra- 
dially stretched by the accreting matter. These field lines 
are located in the disk and above the disk. Most of these 



field lines are connected to the star. The field lines above 
the closed magnetosphere continue to open and reconnect 
in quasi-periodic manner. In case of the lower density 
corona (sec §4), inflation is more efficient, and this leads 
to larger quasi-periodic oscillations of the magnetic flux 
and associated fluxes. 

We analyzed the angular momentum transport between 
the star and the disk and corona. Figure 3b shows the 
distribution of the angular momentum fluxes carried by 
the field fs = rB^'Bp/A'K (left panel) and carried by the 
matter = —prv^Vp (right panel). One can see that for 
r > rA, most of the angular momentum flux is carried by 
the matter (see the region with the high positive angu- 
lar momentum at the right panel). However, for r < 
it is mainly transported by the magnetic field (see left 
panel of Figure 3b). The streamlines in the Figure 3b 
show direction of the flow of angular momentum. One can 
see that matter always carries positive angular momentum 
towards the star, which tends to spin up the star. Mag- 
netic fleld lines carry angular momentum out of the star 
through the field lines threading the area of the funnel 
flow and corona. The situation with the angular momen- 
tum transport seems to be more complicated compared to 
one described by GL79. 

The angular momentum fluxes change with distance 
and with angle. We calculated the angular distribu- 
tion of fluxes Fb (r, 6) = 2nr'^ sin 9iB ■ f and Fm (r, 6) =- 
27rr^ sin Q^m'^ through the spheres of different radii r. Fig- 
ure 4 shows that at a large radius, r = 1.8, matter carries 
most of the angular momentum, while the magnetic field 
also contributes but with the opposite sign. At smaller 
radii, the fluxes become smaller and the largest flux is in 
the region of the funnel flow. At the surface of the star, 
the flux associated with matter is very small. The flux as- 
sociated with the field has two components of the opposite 
sign which cancel each other approximately. 

The distribution of angular momentum flux is more 
complicated compared to that of theoretical models. How- 
ever, the equilibrium state does exist and the ratio rco/TA 
is not very much larger than unity. This means that the 
rotation of the star is efficiently locked to the rotation of 
the inner regions of the disk. 

One can see from Figure 3a, that a signiflcant part of 
the disk is disturbed by the disk-magnetosphere interac- 
tion. Figure 5 shows the distribution of the angular veloc- 
ity of the disk in the equatorial plane. We observed that 
the angular velocity varies around the Keplerian value and 
is ustially slightly smaller then Keplerian. The magnetic 
field is strongly wound up in the disk so that the azimuthal 
field dominates. The inner regions of the disk r < ^ 4 
are appreciably influenced by the magnetic fleld. 

3.3. Dependence on fi and a 

Next, wc took a star in the rotational equilibrium state 
(rco = I-''') and varied the magnetic moment of the star 
fi. Figure 6a shows that when ji increases, the rate of 
change of angular momentum of the star becomes more 
negative, that is a star spins-down. This is connected with 
the fact that for larger /j,. the magnetosphere is larger and 
the flux of the positive angular momentum to the star is 
smaller than in the equilibrium state. For example, for 

= 4, the inner disk radius is at « 1.6. At this in- 
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Fig. 3. — (a). Evolution of the density (color background), and the magnetic field (yellow lines) in rotational equilibrium state for type 
I initial conditions for T = — 200 rotations. The density changes from p = 2 in the disk to p = 0.005 in the corona. The bold red line 
corresponds to /3 = 1 (see §2). (b). Fluxes of angular momentum carried by magnetic field (left panel) and by matter fm (right panel) at 
T = 200. Color background shows value of the fiuxcs, the streamlines with arrows - the direction of the angular momentum flow. 



ner radius, the disk rotates much slower than the star and 
a star strongly spins down. Roughly, the dependence is 




Fig. 4. — Angular distribution of fluxes Fm{r,6) and Fg{r,6) for type I initial conditions. The left-hand panel shows radii along which 
the angular momentum fluxes are calculated. The right-hand panel shows the angular momentum fluxes carried by the magnetic fleld (solid 
lines) and the matter (dashed lines) at different radii r. The numbers show the regions where the angular momentum flux carried by the field 
is carried mostly by the closed (I), radial (II), or open (III) field lines. 
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Lf ~ — 0.02/i + 0.04. The balance between spin-up and 
spin-down torques is similar to that suggest by the GL79 
model, where positive and negative torques are associated 
with regions within and beyond the corotation radius. We 
also observed different regions contributing positive and 
negative angular momentum fluxes. However now the dis- 
tribution of the external to the Alfven surface field lines is 
more complicated. 

We calculated the average value of va for each /i and 
obtained the dependence ~ y^*^, with k = 0.36. This 
coefficient k is somewhat different from that of equation 
(1) where k = 4/7 « 0.57. The difference may be con- 
nected with the fact that in the theoretical analysis the 
magnetic field is assumed to be a pure dipole field every- 
where, whereas the simulations show that the actual field 
is different from a dipole for r > r^- In this and all above 
simulations, the viscosity coefficient was fixed at the value 
a = 0.02. 

We also performed simulation runs for different values 
of the disk viscosity, a = 0, 0.01, 0.02, and 0.03, with 
other parameters fixed. We observed that for larger a 
the accretion rate is larger, and the Alfven radius ta is 
smaller. Thus, the region of positive torque becomes larger 
than the region of negative torque, and i/ is larger. In 
other words, at a higher accretion rate M, incoming mat- 
ter brings more of positive angular momentum, which is 
transferred to larger angular momentum flux carried by 
the field Lf. The angular momentum flux increases with 
a (see Figure 6b) as about Lf ~ a^-^ for a > 0.01. 
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Fig. 6. — (a) The evolution of the angular momentum flux to the 
star carried by the magnetic field L f for different magnetic moments 
of the star fi with other parameters fixed, (b) The evolution of the 
angular momentum flux carried by magnetic field Lf for different 
values of of the disk viscosity a with other parameters fixed. 



4. EQUILIBRIUM STATE OF DISK ACCRETION FOR 
INITIAL CONDITIONS OF TYPE II 





Fig. 5. — Radial distribution of angular velocity of the disk along 
in the equatorial plane at T = 180. The dash line shows Keplerian 
angular velocity. 



200 400 ^ 600 800 
T 

Fig. 7. — The evolution of the matter flux M and the angular 
momentum fluxes carried by matter Lm and magnetic fleld i/ in 
the rotational equilibrium state for the type II initial conditions. 

Here we consider another case corresponding to quite 
different initial conditions referred to as type 11. This case 
has a stronger magnetic field, fi — 10, and a much lower 
coronal density, Pcor = 0.001. The density in the corona 
influences the evolution of the magnetic field in the corona, 
namely, the inflation and opening of the magnetic field 
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Fig. 8. — (a). Evolution of the density (color background), and the magnetic field lines (yellow lines) in rotational equilibrium state of type 
II initial conditions for T = — 250 rotations. The density changes from p = 2 in the disk to p = 0.001 in the corona. The bold red line 
corresponds to /? = 1 (see §2). (b). Fluxes of angular momentum carried by magnetic field fg (left panel) and by matter fm (right panel) at 
T = 250. Color background shows value of the fluxes, the streamlines with arrows - the direction of the angular momentum flow. 



lines is favored. As discussed earlier, the angular momen- turn transport between the disk and star occurs in part 




Fig. 9. — Angular distribution of fluxes Fm{r, 9) and Fg{r, 6) for initial conditions of type II. The left-hand panel shows radii along which 
the angular momenta are calculated. The right-hand panel shows angular momentum fluxes carried by magnetic field (solid lines) and matter 
(dashed lines. The numbers show the regions where the angular momentum flux carried by the field is carried mostly by the closed (I), radial 
(II), or open (III) field lines. 
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through the magnetic field. In addition, we took a dif- 
ferent structure of the disk, which is thinner and located 
initially at = 5, which is closer to the star for such a 
strong field. Wc were able to model the stronger magnetic 
field because wc increased the radius of the star (the in- 
ner boundary) to i?* = 0.5 (versus i?, = 0.35 in other 
simulation runs). In addition, we were able to use a grid 
Nr X Ng = 91x31 which speeded up the simulations and 
allowed longer nms, up to 1000 Pq. 

For this case, we also performed a set of simulations for 
different corotation radii and found the state which gives 
torqueless accretion, L « 0. We found that this state cor- 
responds to a corotation radius Tco ~ 4.6. Figure 7 shows 
the evolution of the fluxes in this case. The angular mo- 
mentum flux associated with the matter is very small as 
in our main case. The flux due to the magnetic field, Lf, 
varies more strongly than in the main case. This is be- 
cause the coronal density is lower so that variations of the 
magnetic field in the corona are larger. In addition, the 
accretion rate M is several times larger than that for our 
main case, because the disk is different compared to the 
main case. 

We analyzed the rotational equilibrium state in detail. 
Figure 8a shows the evolution of the disk and magnetic 
field in the inner part of the simulation region. We ob- 
served that the accretion disk stopped at much larger 
distances, ta ~ 3.1, compared to our main case (r^ ~ 
1.2 — 1.3) because of much stronger magnetic field. This 
radius again coincides with the (3 = 1 line (see the bold 
red line on figure 8a). Wc also estimated the radius ro 
derived from comparison of magnetic and matter torques. 
Again, we obtained similar but a somewhat larger radius, 
ro « 3.2. We observed that in type II case, the ratio 
fco/TA ^ 1-4 — 1.5. This number is slightly larger than 
that in the main case {tco/ta ^ 1-3 — 1.4) which may be 
connected with contribution of the inflated field lines to 
the negative torque. 

Figure 8b shows the angular momentum fluxes carried 
by the field -Fs(r, 9) and by the matter Fm(r, 9) . One can 
see that the disk matter brings positive spin-up torque 
(right panel), while magnetic field carries both negative 
and positive torques (left panel). In this case a significant 
part of the spin-down fiux is carried by the inflated and ra- 
dially stretched field lines. This is different from the type 
I initial conditions case. 

Figure 9 shows angular distribution of ffuxes Fb {r, 6) 
and Fm{r,9) along the spheres of different radii r from 
r = 0.5 to r = 4.5. One can see, that at large distances, 
r = 4.5, there is a large positive angular momentum flux 
carried by the matter of the disk, and there is a large 
negative angular momentum flux carried by the radially 
stretched field lines located above the disk and in the disk. 
Field lines above the disk are inflated field lines which con- 
nect to the disk at large radii r > 20. They carry both, 
positive matter flux, and negative spin-down flux. In this 
region the density is smaller than in the disk, however, 
the field lines are strongly wound and azimuthal velocity 
is large, this is why both fluxes associated with matter 
and the field are large in this region. The plot at r = 3.75 
shows relative input from different sets of field lines. It 
shows that in the region of the closed field lines (region I) 
matter carries positive torque while field lines carry only 



small negative torque. Most of both torque is carried by 
the radially stretched field lines (region II), while the role 
of vertically inflated field lines of "magnetic tower" (region 
III) is very small. At even smaller radii, r = 3,2,1 and 
0.5, both torques become smaller, like in type I case. 

5. EQUILIBRIUM STATE IN APPLICATIONS TO DIFFERENT 

STARS 

The rotational equilibrium state is expected to be the 

most probable state for different disk accreting magnetized 
stars. These include Classical T Tauri Stars, cataclysimic 
variables, and X-ray pulsars. For example, the slowly ro- 
tating CTTSs have spin rates as low as ~ 10% of breakup 
speed (Bouvier et al. 1993). The fact that CTTSs ro- 
tate slowly strongly suggests that they are in rotational 
equilibrium state. 

Below we estimate periods of rotation for different 
weakly magnetized stars. Our simulations were performed 
for cases ta/R* =4 — 6 and can be applied to systems 
with different scales in which this ratio is satisfied. Ob- 
servations show that in CTTS this ratio is in the range 
r^/i?^ = 4-8 (Bouvier et. al., 1993, Edwards et al., 1993), 
or, rA/R* = 3— 10 (Kenyon, Yi & Hartmann, 1996). Sim- 
ilar or smaller magnetospheres are expected in accreting 
millisecond pulsars (van der Klis 2000) and dwarf novae 
cataclysmic variables (Warner 2004). 

For our example we take va/R* = 4, take obtained from 
simulation ratio rco = I-'Ita and derive the angular veloc- 
ity for rotational equilibrium, Oe^: 

3/2 

w 0.09 , (11) 



eg 



K 



R* 



which is about ~ 10% of breakup speed of the star. Here, 
il^K is the Keplerian angular speed at i?*. The value of 
*if we find is close to the one observed in CTTSs. 

5.1. Classical T Tauri Stars (CTTSs) 

The angular velocity of the star can be written as 
ri* = {GM /r^gY^"^ . Combining this with the rotational 
equilibrium condition rco/rA ^ 1-3 — 1.5 (an average 1.4 
was taken here) and using Eqn. 1 for the Alfven radius, 
we can get the critical value of angular velocity Ogg or ro- 
tation period Peq = 2n/^eq- For the CTTS we obtain the 
equilibrium rotation period, 

P.,.4.6d.ys(!^)"'(i2::|^)="x 

where B is the surface magnetic field of the star. Our re- 
sult is in the range of an observational bimodal distribution 
of rotation at period with peaks near 3 and 8 days (At- 
tridge & Herbst 1992). Edwards et al. (1993) proposed 
that the distribution near 8 days may be caused by the 
disk locking. We should note that both peaks may be con- 
nected with the disk locking, but at different parameters, 
say, M and B. 

5.2. Cataclysmic Variables 

We next consider dwarf novae which belong to a sub- 
type of the cataclysmic variables where the magnetic field 
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is expected to be small but still possibly dynamically im- 
portant for the accretion. The accretion disks are stopped 
at small radii by the white dwarf magnetosphere. The ac- 
creting material then leaves the disk and follows the mag- 
netic field lines down to the star's surface in the vicinity 
of the magnetic poles. Taking typical values for the white 
dwarf and accretion disk, we obtain the period of these 
stars in rotational equilibrium state, 



eg 



(■ 



M 
R 



X 18/7 

VlQSG/ V7.0xl08cm7 ' ^^^^ 

The observed periods of Dwarf Novae Oscillations (DNOs) 
are in the range of 7 — 70 s (Warner, 2004). 

5.3. Millisecond Pulsars 

Our simulations are also applicable to millisecond X-ray 
pulsars, which are the accreting low-magnetic-field neu- 
tron stars (van dcr Klis, 2000). The period of such a star 
in the rotational equilibrium is 

-(^)"'(i^)"'^ 

/ B \6/7/ R y8/7 

.lO^cm/ 



IiOSqJ Uo6cmJ ■ ^^^^ 
For all accretion-powered millisecond X-ray pulsars we 

have known, the spin frequencies range from 185Hz to 
435Hz (Wijnands, R. 2004), or, the spin rate is in range 
of 2.3-5.2 ms. These values may be obtained from Eqn.l4 
at smallere M or larger B. 

6. DISCUSSION 

We investigated the conditions for the rotational equi- 
librium or"torqueless" accretion state using axisymmetric 
MHD simulations. In such a state the total angular mo- 
mentum flux to the star is approximately zero. We found 
the equilibrium states by gradually changing the angu- 
lar velocity of the star with the other parameters fixed. 
We considered two main cases: One with relatively low 
magnetic field and dense corona and second with much 
stronger field and lower density corona. We observed that 
in both cases the rotation of the star is approximately 
locked to the rotation of the inner radius of the disk such 



that a star rotates somewhat slower than the inner radius 
of the disk. In the first case the ratio between corota- 
tion radius of the star and the Alfven radius (where 
the disk is disrupted) is rco/rA ^ 1-3 — 1.4. In the second 
case, where stronger inflation of the magnetic field was ob- 
served, tliis ratio is only slightly larger, rco/rA ^ 1.4— 1.5. 
We observed that in the first case, the angular momen- 
tum transport is associated with the closed field lines at 
the inner radii of the disk. Open field lines spin down the 
star, but the role of this spin-down is small. In the sec- 
ond case outflow of angular momentum along the inflated 
field lines is more significant than in the first case, how- 
ever this did not change the result, which probably means 
that the angular momentum transport associated with the 
inner regions of the disk and the region of the funnel flow 
is more significant. Thus, in both cases the magnetic in- 
teraction effectively locks the rotation rate of the star to 
a value which depends mainly on the mass accretion rate 
and the star's magnetic moment. We should note that 
most of coronal region is still matter dominated. This is 
connected with the fact that the magnetic energy-density 
of the dipole decreases with distance as ^ R^^ so that 
it is difficult to set up a magnetically dominated corona. 
Goodson et al. (1997, 1999) were able to model such a 
corona by arranging fast fall-off of coronal density with 
the distance. However, their initial conditions are sufS- 
cicntly far from equilibrium that the torqueless accretion 
was not established. Future simulations with even lower 
coronal densities will help to understand whether a star 
is always locked to the rotation rate of the inner radius of 
the disk. Our simulations for two very different initial con- 
ditions have shown very similar results for the rotational 
equilibrium which may be a sign that the disk-locking may 
be a similar for all slowly rotating stars. 

We applied our simulation results to Classical T Tauri 
stars, where disk locking may explain their slow rotation. 
Also, we estimated the probable periods of rotation of 
other accretion powered systems, such as dwarf novae and 
X-ray millisecond pulsars. 
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